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Abstract: In this study we characterized the relationship between temperature 
and mortality in central Arizona desert cities that have an extremely hot climate. 
Relationships between daily maximum apparent temperature (AT max ) and mortality for 
eight condition-specific causes and all-cause deaths were modeled for all residents and 
separately for males and females ages <65 and >65 during the months May-October for 
years 2000-2008. The most robust relationship was between AT max on day of death and 
mortality from direct exposure to high environmental heat. For this condition-specific 
cause of death, the heat thresholds in all gender and age groups (AT max = 90-97 °F; 32.2- 
36.1 °C) were below local median seasonal temperatures in the study period 
(AT max = 99.5 °F; 37.5 °C). Heat threshold was defined as AT max at which the mortality 
ratio begins an exponential upward trend. Thresholds were identified in younger and older 
females for cardiac disease/stroke mortality (AT max =106 and 108 °F; 41.1 and 42.2 °C) 
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with a one-day lag. Thresholds were also identified for mortality from respiratory diseases 
in older people (AT max = 109 °F; 42.8 °C) and for all-cause mortality in females 
(AT max = 107 °F; 41.7 °C) and males <65 years (AT max = 102 °F; 38.9 °C). Heat-related 
mortality in a region that has already made some adaptations to predictable periods of 
extremely high temperatures suggests that more extensive and targeted heat-adaptation 
plans for climate change are needed in cities worldwide. 

Keywords: apparent temperature; climate; gender; heat-related deaths; hot climate; 
hot cities; temperature threshold 



1. Introduction 

Deaths from exposure to extreme heat have been studied in cities on all inhabited continents [1-7]. 
Motivating this research are high numbers of deaths during extreme heat events and projected 
increases in deaths from heat-related causes due to climate change [8,9]. Heat waves in all types of 
climate zones cause many deaths in cities during short-term atypical weather. In cities with tropical 
and subtropical climates, hot weather also lasts for long stretches during a prolonged warm season 
each year, leading to chronic heat exposure in the population [10,11]. Normal daily temperatures in 
summer are predictably high and hours of exposure to high temperatures in a 24-hour period are also 
longer compared to temperate cities. 

Information about heat-related deaths in hot cities is important to hundreds of millions of people 
who live in them and also contributes to understanding the limits of human tolerance to climatic 
conditions that may become more prevalent worldwide in this century. Average air temperatures in 
most cities are becoming warmer over time for two different reasons: global climate change and the 
urban heat island effect [12,13]. Current temperature-mortality relationships in hot cities are one 
indicator of the extent to which humans have adapted physiologically, behaviorally, and 
technologically to climate and offer insights about living in a future warmer world. 

The purpose of the study was to measure daily mortality response to temperatures during the 
warmest months in southwestern U.S. desert cities that have made adaptations to a very hot climate. 
Our study's focus on enduring heat complements literature on episodic heat events that cause 
extraordinary numbers of deaths in a short time in cities with more temperate climates [3,5-8]. 
We also distinguish among causes of death in hot weather separately by age and gender. The quality of 
mortality data for this type of study is relatively good in the U.S. The world's hottest cities are located 
in deserts of the Middle East and Africa [14], but civil registrations of deaths are low worldwide [15] 
and there are especially high error rates in assigning primary cause of death on death certificates in 
developing nations [16,17]. 
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2. Experimental Section 

2.1. Study Region 

Maricopa County in the Sonoran Desert of south-central Arizona comprises 60% of the state's 
population and it also includes a vast expanse of wilderness (2010 population = 3.8 million, 
area = 9224 sq. mi.). The urbanized area of this county encompasses Phoenix (the largest city), 
15 adjoining municipalities, and three Native American communities (2010 U.S. Census Urban Area 
polygons, total population = 3.6 million, area = 1261 sq. mi.). In the Phoenix area, the 30-year average 
daily maximum temperature during May-October for years 1981-2010 was 99.6 °F, which tied with 
Yuma, Arizona for the highest temperature of all major weather stations in the U.S. [18]. Figure 1 
shows the frequencies and cumulative percentage of days by maximum temperature that occurred 
during our study period. This region also has multi-day summer extreme heat events, which exceed 
historically normal local mean temperatures and usually occur in July and less frequently in June and 
August [19]. 

Figure 1. Frequency and cumulative percentage of days by maximum temperature (°F) at 
the Sky Harbor International Airport Weather Station during the months May-October for 
years 2000-2008. 




Maximum Temperature (°F) 
□ Frequency — Cumulative Percentage | 



Replacement of natural desert and agricultural land cover with residential, commercial, and 
industrial development began in the early 20th century and accelerated during the Sunbelt economic 
expansion from the 1970s to 2008, contributing to a strong urban heat island effect in the settled 
area [20,21]. Ruddell et al. [22] compared historical temperatures between Phoenix (urban) and Gila 
Bend (rural) National Weather Service (NWS) regional surface weather stations in central Arizona and 
found pronounced warming in Phoenix and modest warming in Gila Bend. These differences suggest 
that global climate change is warming natural/rural areas while the combined forces of urbanization 
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and global change are warming urban environments more rapidly. Local and global changes in climate 
have increased the heat burden in the naturally hot desert and are projected to continue doing so in the 
21st century [23]. 

2.2. Data 

The study period covers the warm season in the Arizona desert during the months May-October for 
years 2000-2008. Sky Harbor International Airport in Phoenix is the only NWS first-order surface 
weather station in the study region (first-order stations collect a comprehensive array of weather 
variables every hour on a 24-hour basis). Using hourly air temperature and dew point temperature data 
from Sky Harbor, we calculated hourly apparent temperature (AT) for the study period. We derived 
daily minimum, mean, and maximum AT (AT m ; n , AT mean , AT max ) from the hourly data. Although 
temperatures in urban environments are spatially heterogeneous due to land cover and built 
environment characteristics [24,25], time-series studies often use a single airport or center-city station 
to examine temperature-mortality relationships for a city or metropolitan area (e.g., [4,26]). 

AT is a measure of how hot it feels to people, taking into account relative humidity as well as air 
temperature [27]. We used Equation (1) to estimate hourly AT [28,29]. 

AT = -2.653 + (0.994 x T) + (0.0153 x Tdp 2 ) + wind adj (1) 

where T = Ambient dry-bulb temperature (°C); Tdp = dew point temperature (°C); wind a dj = 
adjustment to AT for given T and wind-speed (meters/second) [30]. 

AT is a good measure for this region because it accounts for human discomfort in the humid 
summer rainy season (July-August) and in the dry warm months. This equation overestimates AT at 
Tdp < 0 °C (32 °F); therefore, for days in our study period with Tdp < 0 °C (<3% of all days), we set 
Tdp = 0 to reduce the error in AT estimates for these days. After this substitution, our AT estimates 
were in close agreement with Steadman's AT values for his temperature -humidity scale (p. 862) [27]. 
In our analysis, temperatures were converted to Fahrenheit (°Celsius x 9/5 + 32 = °Fahrenheit). During 
the study period, median AT max = 99.5 °F; 10% of all days < 83 °F; and 10% of all days > 109 °F. 
Table 1 describes average daily AT max during the study period by month. 



Table 1. Summary statistics for daily maximum apparent temperatures (AT max °F) based 
on daily measurements obtained from the Sky Harbor International Airport Weather 
Station during the months May-October for years 2000 to 2008. 



Month 


AT 

- v 1 max 

Mean 


AT 

- v a max 

SD 


AT 

- v a max 

Minimum 


AT 

■ v 1 max 

Maximum 


May 


92 


9.2 


66 


111 


June 


103 


6.3 


85 


118 


July 


107 


4.6 


93 


118 


August 


105 


5.3 


87 


118 


September 


99 


6.4 


80 


113 


October 


84 


8.3 


60 


104 
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We obtained mortality data for the study period from Part I of death certificates filed in Maricopa 
County. For each record, age, gender, cause of death, and date of death were retrieved from the 
original records. Institutional Review Boards at Arizona State University and the Arizona Department 
of Health Services approved our use of these data for this study. 

Prior literature has estimated heat-related mortality for all causes [1,2,4-6,8,26] and for different 
disease categories, most consistently cardiovascular and respiratory diseases [2,31-34]. To examine 
all-cause mortality related to temperature, we followed the example of others [1,5,26,35,36] and 
excluded most external causes of death. We excluded International Classification of Diseases, 10th 
Revision (ICD-10) codes S00-99, T00-66, T68-98, U00-99, V01-99, W00-99, X00-29, 31, 33-53, 
55-84, Y00-98, Z00-99 and included T67.x, X30, X32 and X54 because they are heat-related. 
We examined eight condition-specific cause-of-death categories derived from a conceptual model of 
medically adverse effects of exposure to high environmental heat and prior literature on environmental 
heat and mortality. These condition-specific cause-of-death categories are described next. The 
Supplementary Material details the reasoning behind selection of categories and Supplemental Table SI 
lists all ICD-10 codes by category: 

Direct exposure to environmental heat. Conditions that might directly cause death in situations of 
high environmental heat were mapped to corresponding primary cause-of-death ICD-10 codes. 
Also included in this category were deaths with terms associated with exposure to high environmental 
heat (e.g., "heat exhaustion") entered as free text in any of the four underlying cause-of-death fields of 
the death certificate Part 1 (Table SI). 

Dehydration. Dehydration is a common and potentially lethal consequence of exposure to high 
environmental heat. The ICD-10 codes for dehydration and its direct consequence — volume 
depletion — and conditions indicative of dehydration (e.g., lack of water) were selected for this 
category. Also included were deaths with "dehydration" entered as free text in any of the cause-of-death 
fields (Table SI). 

Possible consequences of heat and dehydration. This category was defined by ICD-10 codes for 
conditions that are possible consequences of exposure to high environmental heat and/or dehydration 
and terms (e.g., "hypovolemia") entered as free text in the underlying cause-of-death fields (Table SI). 
Figure SI graphically depicts the model that explains our approach to selecting codes considered likely 
to be due to the physiologic and pathophysiologic consequences of heat and dehydration/volume 
depletion/hypovolemia. 

Five other condition-specific categories were selected because prior research has found elevated 
mortality during hot weather. We defined these categories as follows: cardiac disease/stroke (ICD-10 
codes I20.x-I25.x, 149.9, I60.x-I79.x); chronic renal failure (N18.x, N19); heart failure (150. x, 
15 1.6-15 1.9); chronic obstructive pulmonary disease (COPD)Zasthma (J43.x-J45.x, J46,); and other 
respiratory diseases (all other ICD-10 respiratory codes J00-J42, J47-J99.x). 

Mortality at AT max for each cause-specific category and all causes was calculated as the daily ratio 
of actual mortality to baseline May-October ("summer") mortality defined by the equation: 
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Pi t k = = rat i° of number of deaths on year i, month j, day k to average 

' ' Y, (2) 

deaths per day in summer of year i 

where P^ = mortality ratio; Yj^ = number of deaths on year i, month j, day k; Yi - —'^^Y ijk = 

Ni j k 

average deaths per day in summer of year /; N { = number of days in year i from May-October. 
2.3. Statistical Analysis 

We used cubic spline regressions with four knots (equal interval quintiles) to describe the 
temperature-mortality relationship for each category during the study period. The daily time series 
values for AT max , AT mean , and AT min were highly correlated with each other. Consequently, the spline 
regressions for all three AT measures produced very similar results. Models for AT max are presented 
because exposure to peak temperature is most likely during daytime. 

Mortality ratios were the response variables in Equation (3). Predictor variables were AT max with 0, 
1, 2, and 3-day lags as well as an indicator binary variable for heat wave days. Heat waves were 
defined by Meehl and Tebaldi's method [37], which compares each study year to historical 
temperatures and identifies extreme heat events based on three conditions: daily maximum temperature 
must be above the 97.5 percentile of historically normal conditions for at least three days; average 
daily temperature must be above the 97.5 percentile for the entire heat event; and daily maximum 
temperature must be above the 81 percentile for the entire event. The binary variable for heat wave 
days significantly improved model fit for the categories direct exposure to heat, cardio disease/stroke, 
and all-cause deaths. Additional predictor variables for day of the week and public holiday were not 
statistically significant for any category and were eliminated from the models. 

Mortality ratio (Py,*) = spline(AT max . 0 -d ay -iag) + spline(AT max .i. day -ia g ) + 

(3) 

spline(AT max _2-day-i ag ) + spline(AT max . 3 -day-iag) + Pi x heat wave 
Spline function with four knots: 

4 

spline(x) = a 0 + a x x + a 2 x 2 + a 3 x 3 + ^ b m (x - x m ) 3 1(x > x m ) 
x m = the mth knot 

For every condition-specific cause-of-death category and all-cause deaths, we selected the 
temperature metric (AT max with 0, 1, 2 or 3 day lags) that had the largest positive coefficient and a 
statistically significant effect (p < 0.05) on increasing mortality ratio in the spline regression model. 
We then used Poisson regression to model the relationship between the temperature metric and 
mortality in order to identify high temperature thresholds. A lower bound for AT max was established in 
the flat part of the temperature-mortality curve by inspection, and we subsequently tested candidate 
heat thresholds at every integer value of AT max (°F) above the lower bound up to the maximum 
observed temperature in our data (118 °F). From all candidate models, we selected the model and 
associated threshold that provided the best exponential fit as deemed by the goodness-of-fit statistic 
(Pearson chi-squared/df). Our estimate of the heat threshold is the temperature at which the mortality 
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ratio begins an upward exponential trend. We checked our models for overdispersion and found 
significant overdispersion only in the category direct exposure to environmental heat. Hence, for this 
outcome we employed a negative binomial regression model. For each cause-of-death category in 
which a high temperature threshold was identified, we estimated the relative risk (RR) of death above 
the threshold per 1 °F. We note that it was not possible to identify a heat threshold for some of the 
mortality cause-specific categories when the AT max was not a significant predictor of mortality. 

3. Results and Discussion 

3.1. Results 

The statistical analysis in this study is restricted to decedents who were residents of Maricopa 
County. Table 2 displays, for each cause-specific category and all-cause deaths, the number of resident 
deaths (column 2) and the percentage of residents of all deaths that occurred in the county (column 3). 

Figure 2a-i shows the relationships between AT max and mortality ratios (Pij,k in Equation (2)) for 
residents using fitted cubic spline models with 95% confidence intervals (CI) for cause-specific 
categories and all-cause deaths. Corresponding to Figure 2, threshold temperatures and RR/CI for 
deaths above thresholds are reported in Table 2 for category-specific and all-cause deaths with 
statistically significant AT max coefficients (p < 0.05). 

The strongest relationship between mortality and temperature was for direct exposure to high 
environmental heat on day of death (Table 2, row 1 and Figure 2a). In this category, we observed the 
lowest temperature threshold (AT max = 93 °F, lag days = 0) and largest increase in the relative risk of 
death from heat exposure (RR = 1.20); i.e., daily mortality increased by 20% per 1 °F above threshold. 
The relative risk (RR = 1 .05) from possible consequences of heat or dehydration was also increased 
above daily same-day AT max > 106 °F. We found higher temperature thresholds, one-day lags, and 
smaller increases in the relative risk of death from cardiac disease/stroke (AT max = 1 10 °F, RR = 1.03) 
and all-causes (AT max = 107 °F, RR = 1.01) than for deaths directly attributable to heat exposure and 
consequences of heat or dehydration. 

To determine whether heat had different effects on population subgroups, we re-estimated mortality 
models separately for ages <65 years and >65 years (Table 3). We identified significant AT max 
relationships with mortality ratio in both age groups in the categories direct exposure to environmental 
heat and cardiac disease/stroke and in all-cause deaths. For direct exposure to heat, the temperature 
threshold for <65 (AT max = 93 °F) was 3 °F higher than for age >65 (AT max = 90 °F) with zero lag days. 
However, for cardiac disease/stroke, the temperature threshold for age <65 (AT max =106 °F) was 3 °F 
lower than for age >65 (AT max =109 °F) with one-day lags. Similarly, for all-cause deaths, the 
temperature threshold for age <65 (AT max = 1 04 °F) was 3 °F lower than for age >65 (AT max =107 °F) 
with one-day lags. The increases in relative risk above threshold temperatures within each category of 
death were very similar for both age groups. In people age >65, there were statistically significant 
effects of temperature on COPD/asthma deaths (AT max = 1 10 °F, lag days = 0, RR = 1.05). 
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Table 2. Deaths in Maricopa County residents directly caused by exposure to environmental heat and other causes potentially related to high 
temperatures during the months May-October for years 2000-2008. 



Cause of Death 3 


Number of 
County Resident 
Deaths 


Percent Resident 
Deaths of All Deaths c 


Threshold AT max (°F) 
(Lag Days) d 


RR above 
Threshold e 


95%CI f 


Direct exposure to environmental heat 


215 


73.1 


93 (0) 


1.20 


1.16, 1.23 


Dehydration 


409 


93.8 








Possible consequences of heat or dehydration 


1114 


92.1 


106 (0) 


1.05 


1.01, 1.09 


Cardiac disease/stroke 


24,717 


91.2 


110(1) 


1.03 


1.01, 1.05 


Chronic renal failure 


762 


89.8 








Heart failure 


848 


92.0 








Chronic obstructive pulmonary diseases 
(COPD)/asthma 


5342 


92.5 








Other respiratory diseases 


3639 


89.3 








All causes b 


90,284 


91.4 


107(1) 


1.01 


1.00, 1.02 



a See Supplemental Material and Supplemental Table SI for ICD-10 codes and keywords included in cause-of-death categories; b Excludes 9969 deaths of from external 
causes (ICD-10 codes S00-99, T00-66, T68-98, U00-99, V01-99, W00-99, X00-29, 31, 33-53, 55-84, Y00-98, Z00-99) and 187 deaths missing a county of residence 
address; c Number of Maricopa County residents in each cause-of-death category divided by total number of decedents in category; d Coefficient for AT max (appropriate 
lag) has /7-value < 0.05 in best fitting Poisson regression model (Pearson chi squared/df); c Relative risk (RR) is the effect estimate for the increase in mortality per 1 °F 
above the threshold temperature; ' Confidence interval (CI). 
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Figure 2. Temperature-mortality relationships for cause-specific categories and all-cause 
deaths for Maricopa County residents (all ages) during the months May-October for years 
2000-2008. Small black squares ■ are the average mortality ratios (P/^in Equation (2)) or 
ratio of number of actual deaths to average deaths per day aggregated for all years) at a 
given maximum apparent temperature (AT max ). Values < 1.0 indicate fewer than expected 
deaths; values > 1.0 indicate more than expected deaths. Note the scale of mortality ratio 
on the y-axis in (a) is different from the other graphs. Solid black lines are the cubic spline 
regressions predicting daily average mortality ratios at AT max (unsmoothed). Dashed lines 
are the upper and lower 95% confidence intervals (CIs). 
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Figure 2. Cont. 

e) Chronic renal failure f) Heart failure 
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Table 3. Deaths in Maricopa County residents ages <65 years and >65 years directly caused by exposure to environmental heat and other 
causes potentially related to high temperatures during the months May-October for years 2000-2008. 







<65 Years 






>65 Years 




Cause of Death a 


Threshold AT max 
(°F) (Lag Days) c 


RR above 
Threshold d 


95% CI e 


Threshold AT max 
(°F)(Lag Days) c 


RR above 
Threshold" 


95%CI e 


Direct exposure to environmental heat 


93 (0) 


1.19 


1.15, 1.24 


90 (0) 


1.21 


1.15, 1.26 


Possible consequences of heat or 
dehydration 


106 (0) 


1.12 


1.02, 1.23 








Cardiac disease/stroke 


106(1) 


1.03 


1.00, 1.05 


109(1) 


1.03 


1.01, 1.05 


Chronic obstructive pulmonary diseases 
(COPD)/asthma 








110(0) 


1.05 


1.00, 1.10 


All causes b 


104(1) 


1.01 


1.00, 1.02 


107(1) 


1.01 


1.00, 1.02 



a See Supplemental Material and Supplemental Table SI for ICD-10 codes and keywords included in cause-of-death categories; b Excludes deaths from external causes 
(see Table 2); c Coefficient for AT max (appropriate lag) has /?-value < 0.05 in best fitting Poisson regression model (Pearson chi squared/df); d Relative risk (RR) is the 
effect estimate for the increase in mortality per 1 °F above the threshold temperature; e Confidence interval (CI). 



Table 4. Deaths in Maricopa County male residents ages <65 years and >65 years directly caused by exposure to environmental heat and 
other causes potentially related to high temperatures during the months May-October for years 2000-2008. 





Male < 65 Years 




Male > 65 Years 




Cause of Death a 


Threshold AT max 
(°F)(Lag Days) c 


RR above 
Threshold d 


95%CI e 


Threshold AT max 
(°F) (Lag Days) c 


RR above 
Threshold" 


95% CI e 


Direct exposure to environmental heat 


92 (0) 


1.19 


1.15, 1.24 


90 (0) 


1.21 


1.14, 1.29 


Possible consequences of heat or dehydration 


106 (0) 


1.18 


1.06, 1.33 








Other respiratory diseases 








109 (0) 


1.09 


1.01, 1.16 


All causes' 3 


102 (1) 


1.01 


1.00, 1.02 









a See Supplemental Material and Supplemental Table SI for ICD-10 codes and keywords included in cause-of-death categories; b Excludes deaths from external causes 
(see Table 2); c Coefficient for AT max (appropriate lag) has /?-value < 0.05 in best fitting Poisson regression model (Pearson chi squared/df); d Relative risk (RR) is the 
effect estimate for the increase in mortality per 1 °F above the threshold temperature; e Confidence interval (CI). 
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Table 5. Deaths in Maricopa County female residents ages <65 years and >65 years directly caused by exposure to environmental heat and 
other causes potentially related to high temperatures during the months May-October for years 2000-2008. 



Female < 65 Years Female > 65 Years 



Cause-of-Death Category a 


Threshold AT max 
(°F)(Lag Days) c 


RR above 
Threshold d 


95% CI e 


Threshold AT max 
(°F)(Lag Days) c 


RR above 
Threshold" 


95%CI e 


Direct exposure to environmental heat 


97 (0) 


1.19 


1.08, 1.31 


90 (0) 


1.20 


1.12, 1.28 


Possible consequences of heat or dehydration 








105 (0) 


1.05 


1.01, 1.10 


Cardiac disease/stroke 


106(1) 


1.05 


1.01, 1.09 


108 (1) 


1.03 


1.01, 1.05 


Chronic obstructive pulmonary diseases 
(COPD)/asthma 








109 (0) 


1.07 


1.02, 1.13 


All causes b 


107(1) 


1.02 


1.00, 1.04 


107(1) 


1.01 


1.00, 1.02 



a See Supplemental Material and Supplemental Table SI for ICD-10 codes and keywords included in cause-of-death categories; b Excludes deaths from external causes 
(see Table 2); c Coefficient for AT max (appropriate lag) has /?-value < 0.05 in best fitting Poisson regression model (Pearson chi squared/df); d Relative risk (RR) is the 
effect estimate for the increase in mortality per 1 °F above the threshold temperature; e Confidence interval (CI). 
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Re-estimation of separate mortality models for males and females age <65 and >65 showed some 
similarities and differences in heat-related mortality by gender (Tables 4 and 5). All groups had 
relatively low temperature thresholds and high relative risk in the category direct exposure to 
environmental heat (AT max = 90-97 °F, lag days = 0, RR = 1.19-1.21). The relative risk of death from 
the possible consequences of heat or dehydration was also elevated in younger males (AT max =106 °F, 
RR = 1.18) and older females (AT max =105 °F, RR = 1.05) with zero lag days. All groups except 
males >65 had elevated relative risk of death from all causes with one-day lag effects. Table 4 
shows that, above temperature thresholds, the relative risk of death from respiratory diseases 
(non-COPD/asthma) was elevated in males age >65 and Table 5 shows that the relative risk of death 
from COPD/asthma was elevated in females age >65. Relative risk was elevated for cardiac/stroke 
deaths in younger females (AT max = 106 °F, RR = 1.05) and older females (AT max = 108 °F, RR = 1.03) 
with a one-day lag effect. 

3.2. Discussion 

3.2.1. Heat Thresholds and Mortality 

We found increasing mortality above threshold temperatures in Maricopa County residents for 
direct exposure to environmental heat, possible consequences of heat or dehydration, cardiac 
disease/stroke, respiratory illnesses, and all-cause deaths during May-October for years 2000-2008. 
Prior heat and health studies in Maricopa County have demonstrated a significant health response to 
above-normal temperatures, including mortality from heat exposure [38,39] and number of heat-related 
emergency distress dispatches [40]. A national comparative study of heat wave mortality [8] found that 
mortality response was lower in southern cities than in temperate cities but the effect in Phoenix was 
higher than in other southern cities. Yip et al. [38] analyzed mortality risk by age, gender, and race in 
Maricopa County for cause-specific mortality from direct exposure to heat and aggregated 
cardiovascular/respiratory illnesses. Their findings are consistent with ours in important respects as 
described in the next few paragraphs. Our results are also consistent with studies of heat-related 
mortality in many places because deaths associated with high temperatures occurred with very short 
delays [1,2,6,33,35]. Our study contributes new information on protracted threats from heat mortality 
in hot climates and the importance of examining different causes of death in specific population groups. 

In Maricopa County, extremely hot weather is an ever-present danger in summer and the number of 
deaths attributed to heat in Arizona is much higher than in the rest of the U.S. [41]. Although 
heat-mortality threshold temperatures in Arizona desert cities are much higher than in temperate 
cities [2,4,9], people died from direct exposure to environmental heat on days below local median 
seasonal temperatures. In the study period, median daily high AT max was 99.5 °F but the heat tolerance 
threshold was estimated at 93 °F. The NWS advises that prolonged exposure or physical activity 
requires 'extreme caution' when AT (sometimes referred to as the Heat Index) is between 
90 °Fand 105 °F [42] and our estimated daily AT max thresholds for deaths due to direct heat exposure 
are at the low end of that range. 

Prior research suggests that many of the heat exposure deaths in Maricopa County are due to risk 
factors that create acute heat stress situations outdoors, such as working in outdoor occupations [38,43], 
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being homeless [38,43,44], and living in low socioeconomic status and hotter neighborhoods [45]. 
Lack of air conditioning or non-operational units on hot days also contribute to indoor heat 
deaths [38,46]. Temperature effects on direct heat exposure deaths were statistically significant for all 
age and gender groups in Maricopa County: younger and older males and females had similarly 
elevated relative risks above relatively low temperature thresholds in this climate throughout the 
summer months. 

Our results for direct heat exposure deaths were robust and because of their irrefutable connection 
to heat, they bear closer examination. We identified 215 deaths in residents (average 25.5 
deaths/season) in nine summers (2000-2008) using all available information (ICD-10 code and 
keywords) in death certificates Part I. Mortality in males age <65 outnumbered all other groups: males 
<65 = 112 deaths (52.1%); females <65 = 22 (10.2%); males >65 = 45 (20.9%); females >65 = 36 
(16.7%). Yip et al. [38] reported that the majority of deaths attributed directly to heat during the 2005 
heat wave were in white males with median age 55, and the majority of those deaths occurred outdoors. 
During five summers (2007-2011) of operating a heat surveillance system that uses multiple sources 
of information to investigate deaths suspected of being associated with environmental heat, the 
Maricopa County Department of Public Health confirmed 278 heat-associated deaths in residents 
(average 55.6 deaths/season) [39]. In all likelihood, our study undercounted the number of deaths 
associated with direct exposure to environmental heat. Studies that use only the cause-of-death ICD-10 
code count even fewer cases. 

Most literature heavily emphasizes sensitivity of the elderly to environmental heat, especially from 
cardiovascular and respiratory causes [1,28,32,34,36,47,48]. Heat effects on male and female mortality 
appear to vary by cause of death and study location [1]. Some studies have not found any gender 
differences [47,48] but others found the effect of heat mortality was stronger for females [5,36,49,50]. 
It is important to distinguish among different causes of death in order to understand differences in 
heat-related mortality between the young and old or females and males. 

The only consistent age difference in our results was elevated relative risk above the temperature 
thresholds from respiratory illnesses for older males and females. Unlike most prior studies, with the 
exception of respiratory illnesses, we did not find consistent differences in the heat-related relative 
risks for other categories of death by age. If fact, our analysis showed that threshold temperatures were 
lower in people <65 years than in the elderly for cardiac disease/stroke and all-causes deaths. Mortality 
caused by possible consequences of heat or dehydration had significant relationships with daily AT max 
for younger males and older females. 

The relative risk of death from cardiac disease/stroke above the heat threshold was elevated in 
younger and older females but not at all in males of any age. We found a significantly different 
temperature-mortality relationship between males and females for cardiovascular disease, which could 
suggest gender-specific physiological differences in the effects of sustained hot weather on heat stroke 
in men and women (e.g., [51]). 

Cardiac disease/stroke deaths are numerically consequential because they alone accounted for 
27.4% of all-cause deaths in county residents. The NWS guidelines indicate that heat-related illnesses 
are likely when AT is in the 'danger' zone between 105 and 130 °F [42], and our estimated heat 
thresholds for death from cardiac disease/stroke, respiratory diseases and other possible consequences 
of heat and dehydration fall at the lower end of this range between AT max =105 and 1 10 °F. 
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We observed a one-day lag between daily high temperature and cardiac disease/stroke deaths and 
all-cause deaths. Yip et al. [38] also found a significant association of temperature with deaths in a 
combined cardiovascular/respiratory category. Using time series analysis for the months June-September 
in 2000-2005 (and focusing especially on a heat wave period in 2005), those authors found an 
exponential Heat Index-mortality relationship with estimated RR = 1.06 (95% CI: 1.00, 1.13) for each 
1 °F increase in temperature. This is in the middle of our estimates for cardiac/stroke, COPD/asthma, 
and other respiratory diseases (RR = 1.03-1.09) in Tables 2 through 5. Yip et al. [38] did not estimate 
thresholds but they reported that cardiovascular/respiratory mortality increased exponentially with 
temperature as the Heat Index "approaches 100 °F." Based on the evidence, we conclude that 
prolonged high temperature during the summer in hot cities is very likely related to higher mortality 
from cardiovascular disease. 

Chronic heat may increase the relative risk of death for people in hot climates because the body's 
thermoregulatory response to heat stress includes increased skin blood flow and peripheral blood 
pooling, which place a burden on the heart to maintain adequate cardiac output and blood pressure. 
Increased sweating can lead to dehydration and hemoconcentration, which then causes coronary 
thrombosis and stroke. Our approach to categorizing cardiovascular deaths differs in some respects 
from the approach taken by others, who have generally done their analyses based on grouping causes 
of death according to blocks of ICD-10 codes without considering the underlying pathophysiologic 
mechanism that would link exposure to high environmental temperature and the condition listed as the 
cause of death. 

We categorized some deaths that are in the "diseases of the circulatory system" block of codes 
(100-199; e.g., deaths due to arrhythmias and sudden death) with deaths in the "endocrine, nutritional 
and metabolic diseases" block of codes (E00-E90; e.g., hyperosmolality and hypernatremia) because 
the effect of heat in causing a death attributed to the condition is both direct and immediate and the 
effect of heat might be expected to be large for these conditions. We did find a lower threshold 
temperature and a shorter lag for mortality in this "possible consequences of heat and dehydration" 
category than in our "cardiac disease/stroke" category (Table 2). We analyzed deaths due to heart 
failure (150), which are included in the "diseases of the circulatory system" block in most analyses, as 
a separate cause of death because the additional cardiac burden that accompanies heat stress might be 
expected to be especially dangerous for this condition. We did not find heart failure mortality to be 
significantly associated with mortality risk in this setting. 

Several methodologies using time-series techniques have been employed in prior studies to 
characterize temperature-mortality relationships [2,4,52]. For example, Curriero et al. [2] analyzed the 
temperature-mortality spectrum using relative risk functions for 1 1 US cities (which did not include 
Phoenix) during the period 1973-1994. Curriero et a/.'s analysis yielded approximately J-shaped 
relative-risk curves and estimated the so-called minimum mortality temperature (MMT) that separates 
mortality associated with the "cold slope" (temperatures lower than the MMT) from the "hot slope" 
mortality (temperatures higher than the MMT). Similarly, Muggeo (52) developed a unified 
model-based approach to identify break points in mortality data as a function of temperature that could 
indicate increasing physical risk from cold and heat. Muggeo 's approach was used in Milan, Italy, 
which is characterized by a similar J-shape to those identified by Curriero et al. Our statistical method 
for estimating heat thresholds was adapted from McMichael et a/.'s [4] TSOTHURM' study that 
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quantified mortality relationships with cold and heat in a dozen world cities. Due to the absence of a 
"cold season" in Maricopa County, we focused only on warm months to reduce seasonal fluctuations 
in mortality. This approach is well suited for settings in hot climates, such as the desert southwestern 
U.S., where the temperature -mortality relationship includes a substantial temperature interval with a 
"flat" mortality ratio before a significant mortality increase is observed. It is worth noting that our heat 
threshold estimates provide the value of temperature at which mortality ratio begins an exponential 
trend, which should not be confused with the turning points derived in some previous studies. 

In all retrospective studies of temperature-associated mortality, the magnitudes of thresholds for 
categories of deaths indirectly related to heat must be interpreted with caution. The ISOTHURM study 
states that threshold estimates can be imprecise when "data are sparse over a very narrow temperature 
range" [4]. A time series of mortality records longer than nine years would increase the number of 
observed days that deviate widely from local normal temperatures (Figure 1) and reduce uncertainties 
around mortality estimates. Figure 2b— i show high variability in mortality on unseasonably cool days 
(AT max < 78° F). Moreover, high variability in mortality was observed in several death categories (e.g., 
Figure 2e-h) on extremely hot days. This may not be surprising because relatively few extreme days 
(e.g., AT max > 1 14 °F, n = 22) occurred during our study period. 

Estimating thresholds for categories of deaths that are not directly attributable to heat is subject to 
measurement error because a heat effect can be masked by the majority of deaths in that category that 
are not heat-related. The temperature threshold estimates are much higher for categories of deaths 
indirectly associated with heat than for deaths directly caused by heat and this may be partly due to 
overestimating the thresholds for indirect causes. The small numbers of deaths in some cause-of-death 
categories and population subgroups also warrant caution in interpreting the significance of 
temperature thresholds and in comparing differences between age and gender groups. There were 
relatively few deaths in the categories for which we did not identify temperature thresholds (with the 
exception of direct heat exposure). Consequently, some differences in observed relationships between 
population subgroups and cause-specific categories of deaths are affected by significant uncertainty 
associated with reduced statistical power for estimates that were based on small death counts. While 
the threshold measured for each group may be subject to some measurement error, we did find several 
temperature effects on deaths in younger people and this should be enough to justify more thorough 
examinations of the presumption that elderly people are at greater risk of heat-related deaths in 
hot climates. 

Other limitations of this study include the following. Air quality was not measured in our study so 
we do not know if there are interactions between heat and air pollution. Yip et al. [38] found weak 
associations of heat deaths with PMio and ozone in Maricopa County. As noted in Section 2.2, 
exposure misclassification is potentially a problem in studies of temperature-related mortality in large 
urban areas because it is common to use one weather station or to average temperatures over a handful 
of regional stations. Meteorological measurements at one place may be different than measurements in 
others at a given time due to the heat island effect, differences in elevation, or wind circulation. The 
regional temperature-mortality threshold estimates in this study are based on meteorological conditions 
at the Sky Harbor International Airport station located in the central city. Air temperatures are higher 
there compared to suburban neighborhoods but not higher than some inner-city neighborhoods and 
urban fringe areas during the hottest part of the day [24]. There are also different levels of 
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socioeconomic vulnerability to heat stress in populations and neighborhoods throughout this urbanized 
area [53]. Thus, temperature thresholds for heat-related deaths may be lower in some places and higher 
in others due to intra-urban differences in physical environments and population characteristics. 
A recent study used spatial methods to interpolate maximum daily temperature measurements from 
19 weather stations in Brisbane, Australia and found that time series models using data from the 
center-city station and a spatiotemporal model of distributed temperatures performed equally well and 
provided similar estimates of short-term heat mortality for the whole region [54]. It is beyond the 
scope of our study to investigate geographic differences in thresholds. 

3.2.2. Implications of Climate Change for Future Heat Mortality 

Climate projections for central Arizona indicate that future air temperatures will continue to rise in 
the 21st century due to the effects of global climate change and continuing urban expansion. The North 
American Regional Climate Change Assessment Program (NARCCAP) provides an ensemble of 
climate change simulations that couple several regional climate models with atmospheric-ocean 
general circulation models that are driven by the IPCC SRES A2 scenario [54,56]. Analyzing all 
NARCCAP model combinations for central Arizona, Grossman-Clarke et al. [57] concluded that 
regional mean daily maximum air temperature in mid-century (averaged for all 2041-2070 summer 
seasons) is projected to be 110.1 °F. NARCCAP simulations and the Grossman-Clarke et al. analysis 
did not consider the effects of potential urban expansion on future temperature. During the past several 
decades, however, evening and nighttime temperatures in the city of Phoenix increased relative to rural 
sites [22]. Georgescu et al. [23,58] have demonstrated that projected expansion of urban land uses will 
induce higher mid-century summer temperatures in the region. Thus, higher summer temperatures are 
expected in the future because of both global and local forcings, and this could dramatically amplify 
future heat-related mortality in Arizona. 

In the past, humans have adapted physiologically, behaviorally, and technologically to a broad 
range of temperatures through measures such as better health care and housing, weather warnings, and 
emergency response [59]. The major technological adaptation to hot weather in the U.S. is the 
invention and adoption of air conditioning, which accounted for 80% of the substantial national decline 
in heat-related deaths during the 20th century [60]. Urban mortality studies have shown that air 
conditioning is an important preventative factor against heat deaths [61]. It is likely that air 
conditioning has forestalled higher heat-related mortality in central Arizona because its widespread 
adoption preceded this region's rapid urbanization in the latter half of the 20th century. Today, 97% of 
all occupied housing units in the Phoenix metropolitan area have central air compared to 66% of 
occupied units in the nation [62]. Coverage is 92% in households below the poverty line [62], although 
the costs of electricity and repairing broken units are barriers to operating air conditioners in many 
low-income Phoenix households [63]. Control of indoor temperatures appears to have nearly reached 
the saturation point in this region, leaving accessibility gaps only in low-income households. However, 
hotter summer weather and continuing drought in the Southwest [64] could increase the potential for 
regional power failures by increasing demand for electricity and interrupting the water supply to 
power plants. 
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Cities in all types of climate regimes are becoming warmer in the 21st century [12,13] and may be 
subjected to more Arizona-like weather without the benefit of widespread air conditioning that 
Arizonans use to cope with extreme heat. A recent study asked whether humans will be able to adapt to 
the future climate in the event that projected levels of global warming would exceed the physiological 
limit to intolerably high heat stress [65]. Sherwood and Huber [65] maintain that global warming will 
make some populous areas of the world uninhabitable if the average global temperature rises by at 
least 7 °C (12.6 °F). Even if future temperatures do not exceed the limit of human physiological 
tolerance, numerous analyses predict that global heat-related deaths in cities are very likely to increase 
in the future. High heat mortality in central Arizona, a region that has already made some adaptations 
to predictable periods of extremely high temperatures, demonstrates that more extensive heat-adaptation 
plans for climate change are needed in cities worldwide. 

4. Conclusions 

This study found significant temperature-mortality relationships in cities with an extremely hot 
climate during the warmest season when median daily AT max = 99.5 °F. The temperature thresholds, or 
AT max at which the mortality ratio begins an upward exponential trend varied between a low of 
AT max = 93 °F (zero-day lag) for mortality due to direct exposure to environmental heat to a high of 
AT max = 110 °F (one-day lag) for cardiac disease/stroke. All age and gender subgroups had elevated 
relative risk above a threshold for direct heat exposure deaths. Females of all ages had elevated relative 
risk for deaths from cardiac disease/stroke above threshold temperatures. On the hottest days, relative 
risk was elevated for deaths from respiratory diseases in the elderly and for all-cause deaths in females 
and younger males. These results suggest it is important to look beyond analyses of all-cause deaths 
and extreme heat events in order to understand the different consequences of prolonged heat stress for 
vulnerable populations. More detailed analyses, such as this one, are required to strengthen 
multi-faceted adaptive measures that address all kinds of heat situations. In view of projected future 
climate change from global warming and urban expansion, it appears likely that there will be amplified 
heat-health risks in cities worldwide and managing risks under future conditions will require 
adaptation planning for chronically high temperatures. 

Acknowledgements 

This research was supported by grants from the National Science Foundation (GEO-08 16168, 
BCS-1026865, and EF-1049251). The Arizona State Department of Health Services provided mortality 
data for the study. The authors gratefully acknowledge the assistance of Vjollca Berisha and Sarah 
Santana, Office of Epidemiology, Maricopa County Department of Public Health. We are indebted to 
our multi-disciplinary research team on the Urban Vulnerability to Climate Change project, especially 
Susanne Grossman-Clarke, and to David Hondula and Matei Georgescu for their advice. We benefited 
greatly from the insightful comments of three anonymous reviewers, which helped us to improve the 
manuscript. 



Int. J. Environ. Res. Public Health 2014, 11 



3322 



Author Contributions 

Sharon Harlan, Gerardo Chowell, and Diana Petitti designed the study and supervised the coding 
and analysis of data. Harlan conducted the literature review and wrote the paper. Chowell designed the 
study's analytic strategy. Shuo Yang analyzed the data and prepared tables and figures. Petitti designed 
the classification system for selection of conditions and diseases related to exposure to high 
environmental temperatures. Emmanuel Morales Butler programmed the disease classification system 
and contributed to supplementary material. Benjamin Ruddell and Darren Ruddell contributed to the 
intellectual content of the paper. Darren Ruddell also obtained and analyzed temperature data. 
All authors reviewed and commented on the manuscript. 

Conflicts of Interest 

The authors declare no conflict of interest. 
References 

1. Bell, M.L.; O'Neill, M.S.; Ranjit, N.; Borja-Aburto, V.H.; Cifuentes, LA.; Gouveia, N.C. 
Vulnerability to heat-related mortality in Latin America: A case-crossover study in Sao Paulo, 
Brazil, Santiago, Chile and Mexico City, Mexico. Int. J. Epidemiol. 2008, 37, 796-780. 

2. Curriero, F.C.; Heiner, K.S.; Samet, J.M.; Zeger, S.L.; Strug, L.; Patz, JA. Temperature and 
mortality in 1 1 cities of the eastern United States. Am. J. Epidemiol. 2002, 155, 80-87. 

3. D'Ippoliti, D.; Michelozzi, P.; Marino, C; de'Donato, F.; Menne, B.; Katsouyanni, K.; 
Kirchmayer, U.; Analitis, A.; Medina-Ramon, M.; Paldy A.; et al. The impact of heat waves on 
mortality in 9 European cities: Results from the EuroHEAT project. Environ. Health 2010, 9, 37, 
doi:10.1186/1476-069X-9-37. 

4. McMichael, A.J.; Wilkinson, P.; Kovats, R.S.; Pattenden, S.; Hajat, S.; Armstrong, B.; 
Vajanapoom, N.; Niciu, E.M.; Mahomed, H.; Kingkeow, C; et al. International study of 
temperature, heat and urban mortality: The "ISOTHURM" project. Int. J. Epidemiol. 2008, 37, 
1121-1131. 

5. Son, J.Y.; Lee, J.T.; Anderson, G.B.; Bell, M.L. The impact of heat waves on mortality in seven 
major cities in Korea. Environ. Health. Perspect. 2012, 120, 566-571. 

6. Tong, S.; Wang, X.Y.; Barnett, G.A. Assessment of heat-related health impacts in 
Brisbane, Australia: Comparison of different heat wave definitions. PLOS ONE 2010, 5, 
doi:10.1371/journal.pone.0012155. 

7. Egondi, T.; Kyobutungi, C; Kovats, S.; Muindi, K; Ettarh, R.; Rocklov, J. Time-series analyses 
of weather and mortality patterns in Nairobi's informal settlements. Global Health Action 2012, 5. 
Available online: http://www.globalhealthaction.net/index.php/gha/article/view/19065 (accessed 
on 23 January 2014). 

8. Anderson, G.B.; Bell, M.L. Heat waves in the United States: Mortality risk during heat waves and 
effect modification by heat wave characteristics in 43 U.S. communities. Environ. Health Persp. 
2011, 779,210-218. 



Int. J. Environ. Res. Public Health 2014, 11 



3323 



9. Hajat, S.; Kosatsky, T. Heat-related mortality: A review and exploration of heterogeneity. 
J. Epidemiol. Commun. H. 2010, 64, 753-760. 

10. Peretz, C; Biggeri, A.; Alpert, P.; Baccini, M. The Effect of Heat Stress on Daily Mortality in Tel 
Aviv, Israel. In National Security and Human Health Implications of Climate Change; Fernando H.J., 
Klaic Z.B., McCulley J.L., Eds.; Springer: New York, NY, USA, 2012; pp. 241-251. 

11. Balogun, A.A.; Balogun, LA.; Adeyewa, Z.D. Comparisons of urban and rural heat 
stress conditions in a hot-humid tropical city. Glob. Health Action 2010, 3, 5614, 
doi:10.3402/gha.v3i0.5614. 

12. Contribution of Working Group I to the Fourth Assessment Report of the Intergovernmental Panel 
on Climate Change; Solomon, S., Qin, D., Manning, M., Chen, Z., Marquis, M., Averyt, K.B., 
Tignor, M., Miller, H.L., Eds.; Cambridge University Press: Cambridge, UK, 2007. 

13. Grimmond, S.B. Urbanization and global environmental change: Local effects of urban warming. 
The Geogr. J. 2007, 173, 83-88. 

14. Osborn, L. Hottest Cities in the World. Current Results: Research News & Science Facts. 
Available online: http://www.currentresults.eom/W eather-Extremes/hottest-cities-in-the-world.php 
(accessed on 16 December 2013). 

15. World Health Organization. Civil Registration: Why Counting Births and Deaths is Important. 
Fact Sheet No. 324. April 2013. Available online: http://www.who.int/mediacentre/ 
factsheets/fs324/en/ (accessed on 16 December 2013). 

16. Haque, A.S.; Shamim, K.; Siddiqui, N.H.; Irfan, M.; Khan, JA. Death certificate completion 
skills of hospital physicians in a developing country. BMC Health Serv. Res. 2013, 13, 
doi:10.1 186/1472-6963-13-205. 

17. Rao, C; Bradshaw, D.; Mathers, CD. Improving death registrations and statistics in developing 
countries: Lessons from sub-Saharan Africa. S. Afr. J. Demogr. 2004, 9, 81-99. 

18. National Oceanic and Atmospheric Administration Comparative Climatic Data for the United 
States through 2012. Available online: http://wwwl.ncdc.noaa.gov/pub/data/ccd-data/CCD-2012.pdf 
(accessed on 8 January 2014). 

19. Grossman-Clarke, S.; Zehnder, J.A.; Loridan, T.; Grimmond, S.B. Contribution of land use 
changes to near-surface air temperatures during recent summer extreme heat events in the Phoenix 
metropolitan area. J. Clim. 2010, 49, 1649-1664. 

20. Chow, W.T.L.; Brennan, D.; Brazel, A.J. Urban heat island research in Phoenix, Arizona: 
Theoretical contributions and policy applications. Bull. Amer. Meteor. Soc. 2012, 93, 517-530. 

21. Brazel, A.; Selover, N.; Vose, R.; Heisler, G. The tale of two climates — Baltimore and Phoenix 
urban LTER sites. Clim. Res. 2000, 15, 123-135. 

22. Ruddell, D.; Hoffman, D; Ahmad, O.; Brazel, A. Historical threshold temperatures for Phoenix 
(Urban) and Gila Bend (Desert), Central Arizona, USA. Clim. Res. 2013, 55, 201-215. 

23. Georgescu, M.; Moustaoui, M.; Mahalov, A.; Dudhia, J. Summer-time climate impacts of 
projected megapolitan expansion in Arizona. Nat. Clim. Chang. 2012, 3, 37-41. 

24. Ruddell, D.M.; Harlan, S.L.; Grossman-Clarke, S.; Buyantuyev, A. Risk and Exposure to Extreme 
Heat in Microclimates of Phoenix, AZ. In Geospatial Techniques in Urban Hazard and Disaster 
Analysis; Showalter, P., Lu, Y., Eds; Springer: New York, NY, USA, 2010; pp. 179-202. 



Int. J. Environ. Res. Public Health 2014, 11 



3324 



25. Zhang, K.; Oswald, E.M.; Brown, D.G.; Brines, S.J.; Gronlund, C.J.; White -Newsome, J.L.; 
Rood, R.B.; O'Neill, M.S. Geospatial exploration of spatial variation of summertime temperatures 
in the Detroit metropolitan region. Environ. Res. 2011, 111, 1046-1053. 

26. Braga, A.L.F.; Zanobetti, A.; Schwartz, J. The time course of weather-related deaths. 
Epidemiology 2001, 12, 662-667. 

27. Steadman, R.G The assessment of sultriness. Part I: A temperature-humidity index based on 
human physiology and clothing science. J. Appl. Meteor. 1979, 18, 861-873. 

28. Basu, R. High ambient temperature and mortality: A review of epidemiologic studies from 
2001-2008. Environ. Health 2009, 8, 1-13. 

29. Kalkstein, L.S.; Valimont, K.M. An evaluation of summer discomfort in the United States using a 
relative climatological index. Bull. Amer. Meteor. Soc. 1986, 67, 842-848. 

30. Steadman, R.G. The assessment of sultriness. Part II: Effects of wind, extra radiation and 
barometric pressure on apparent temperature. J. Appl. Meteor. 1979, 18, 874-885. 

31. Barnett, A.G Temperature and cardiovascular deaths in the U.S. elderly: Changes over time. 
Epidemiology 2007, 18, 369-372. 

32. Basu, R.; Dominici, F.; Samet, J.M. Temperature and mortality among the elderly in the United 
States: A comparison of epidemiologic methods. Epidemiology 2005, 16, 58-66. 

33. Braga, A.L.F.; Zanobetti, A.; Schwartz, J. The effect of weather on respiratory and cardiovascular 
deaths in 12 U.S. cities. Environ. Health Perspect. 2002, 110, 859-863. 

34. Vaneckova, P.; Beggs, P.J.; de Dear, R.J.; McCracken, K.W. Effect of temperature on mortality 
during the six warmer months in Sydney, Australia between 1993 and 2004. Environ. Res. 2008, 
108, 361-369. 

35. Basu, R.; Feng, W.; Ostro, B.D. Characterizing temperature and mortality in nine California 
counties, 1999-2003. Epidemiology 2008, 19, 138-145. 

36. Medina-Ramon, M.; Zanobetti, A.; Cavanaugh, D.P.; Schwartz, J. Extreme temperatures and 
mortality: Assessing effect modification by personal characteristics and specific cause of death in 
a multi-city case-only analysis. Environ. Health Perspect. 2006, 114, 1331-1336. 

37. Meehl, G.A.; Tebaldi, C. More intense, more frequent, and longer lasting heat waves in the 21st 
century. Science 2004, 305, 994-997. 

38. Yip, F.Y.; Flanders, W.; Wolkin, A; Engelthaler, D.; Humble, W.; Neri, A.; Lewis, L.; Backer, L.; 
Rubin, C. The impact of excess heat events in Maricopa County, Arizona: 2000-2005. Int. J. 
Biometeorol. 2008, 52, 765-772. 

39. Maricopa County Department of Public Health. Annual Heat Reports 2007-2011. Available 
online: http://www.maricopa.gov/publichealth/Services/EPI/Reports/heat.aspx (accessed on 17 
October 2013). 

40. Golden, J.; Hart, D.; Brazel, A.; Luber, G; Phelan, P. A biometeorology study of climate and 
heat-related morbidity in Phoenix from 2001 to 2006. Int. J. Biometeorol. 2008, 52, 471-480. 

41. Centers for Disease Control and Prevention (CDC). Heat-related deaths — United States, 
1999-2003. Morb. Mortal Wkly. Rep. 2006, 55, 796-798. 

42. National Oceanic and Atmospheric Administration (NOAA). Heat Wave: A Major Summer Killer. 
U.S. Department of Commerce. Available online: http://www.nws.noaa.gov/om/brochures/ 
heatwave.pdf (accessed on 17 October 2013). 



Int. J. Environ. Res. Public Health 2014, 11 



3325 



43. Petitti, D.B.; Harlan, S.L.; Chowell-Puente, G.; Ruddell, D. Occupation and environmental 
heat-associated deaths in Maricopa County, Arizona: A case-control study. PLOS ONE 2013, 8, 
e62596, doi:10.1371/journal.pone.0062596. 

44. Sanchez, C. Tricks of the Shade: Heat-related Coping Strategies of Urban Homeless Persons in 
Phoenix, Arizona. Master Thesis. Arizona State University, Tempe, AZ, USA, 201 1. 

45. Harlan, S.L.; Brazel, A.; Prashad, L.; Stefanov, W.L.; Larsen, L. Neighborhood microclimates and 
vulnerability to heat stress. Soc. Sci. Med. 2006, 63, 2847-2863. 

46. Maricopa County Department of Public Health. Annual Heat Reports 2012. Available online: 
http://www.maricopa.gov/publichealth/Services/EPI/pdf/heat/2012annualreport.pdf (accessed on 
17 October 2013). 

47. Basu, R.; Ostro, B.D. A multicounty analysis identifying the populations vulnerable to mortality 
associated with high ambient temperature in California. Am. J. Epidemiol. 2008, 168, 632-637 '. 

48. Choi, G.; Choi, J.; Kwon, H. The impact of high apparent temperature on the increase of 
summertime disease-related mortality in Seoul: 1991-2000. J. Prev. Med. Public Health 2005, 38, 
283-290. 

49. Diaz, J.; Jordan, A.; Garcia, R.; Lopez, C; Alberdi, J.C.; Hernandez, E.; Otero, A. Heat waves in 
Madrid 1986-1997: Effects on the health of the elderly. Int. Arch. Environ. Occup. Health 2002, 
75, 163-170. 

50. Rainham, D.G.C.; Smoyer-Tomic, K.E. The role of air pollution in the relationship between a heat 
stress index and human mortality in Toronto. Environ. Res. 2003, 93, 9-19. 

51. Grogan, H.; Hopkins, P.M. Heat stroke: Implications for critical care and anaesthesia. Br. J. 
Anaesth. 2002, 88, 700-707. 

52. Muggeo, V.M.R. Modeling temperature effects on mortality: Multiple segmented relationships 
with common break points. Biostatistics 2008, 9, 613-620. 

53. Harlan, S.L.; Declet-Barreto, J.H.; Stefanov, W.L.; Petitti, D.B. Neighborhood effects on heat 
deaths: Social and environmental predictors of vulnerability in Maricopa County, Arizona. 
Environ. Health Persp. 2013, 121, 197-204. 

54. Guo, Y.; Barnett, A.G.; Tong, S. Spatiotemporal model or time series model for assessing 
city-wide temperature effects on mortality? Environ. Res. 2013, 120, 55-62. 

55. Mearns, L.O.; Arritt, R.; Biner, S.; Bukovsky, M.S.; McGinnis, S.; Sain, S.; Caya, D.; Correia, J., Jr.; 
Flory, D.; Gutowski, W.; et al. The North American Regional Climate Change Assessment 
Program: Overview of Phase I results. BAMS 2012, 93, 1337-1362. 

56. North American Regional Climate Change Assessment Program (NARCCAP). Available online: 
http://www.narccap.ucar.edu/about/index.html (accessed on 18 December 2013). 

57. Clarke, S.G.; Clarke, T.R.; Harlan, S.L. Extreme summer heat in Phoenix, Arizona (USA) under 
global climate change (2140-2070). Under review. 

58. Georgescu, M.; Moustaoui, M.; Mahalov, M.; Dudhia, J. An alternative explanation of the 
semiarid urban area "oasis effect". J. Geophys. Res. Atmos. 2011, 116, D24113. Available online: 
http://dx.doi.org/10.1029/2011JD016720 (accessed on 23 January 2014). 

59. Davis, R.E.; Knappenberger, P.C.; Novicoff, W.M.; Michaels, P.J. Decadal changes in human 
heat-related mortality in the eastern United States. Clim. Res. 2002, 22, 175-184. 



Int. J. Environ. Res. Public Health 2014, 11 



3326 



60. Barreca, A.; Clay, K.; Deschenes, O.; Greenstone, M.; Shapiro, J.S. Adapting to Climate Change: 
The Remarkable Decline in the U.S. Temperature-mortality Relationship over the 20th Century. 
Working Paper 18692. National Bureau of Economic Research. 2013. Available online: 
http://www.nber.org/papers/wl8692 (accessed on 17 October 13). 

61. Bouchama, A.; Dehbi, M.; Mohamed, G.; Matthies, F.; Shoukri, M.; Menne, B. Prognostic factors 
in heat-wave related deaths: A meta-analysis. Arch. Intern. Med. 2007, 167, 2170-2176. 

62. U.S. Bureau of the Census. 2011 American Community Housing Survey Metropolitan Areas. 
Table C03AOM. Available online: http://www.census.gov/housing/ahs/data/2011ma.html 
(accessed on 17 October 2013). 

63. Hayden, M.H.; Brenkart- Smith, H.; Wilhelmi, O.V. Differential adaptive capacity to extreme heat: 
A Phoenix, Arizona case study. Weather Clim. Soc. 2011, 3, 269-280. 

64. Overpeck, J.; Udall, B. Dry times ahead. Science 2010, 328, 1642-1643. 

65. Sherwood, S.C.; Huber, M. An adaptability limit to climate change due to heat stress. P. Natl. 
Acad. Sci. USA 2010, 107, 9483-9484. 

© 2014 by the authors; licensee MDPI, Basel, Switzerland. This article is an open access article 
distributed under the terms and conditions of the Creative Commons Attribution license 
(http://creativecommons.Org/licenses/by/3.0/). 



